Quantum dot ground state energies and spin polarizations: soft versus hard chaos 
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We consider how the nature of the dynamics affects ground state properties of ballistic quantum dots. We 
find that "mesoscopic Stoner fluctuations", that arise from the residual screened Coulomb interaction, are very 
sensitive to the degree of chaos. It leads to ground state energies and spin-polarizations whose fluctuations 
strongly increase as a system becomes less chaotic. The crucial features are illustrated with a model that depends 
on a parameter that tunes the dynamics from nearly integrable to mostly chaotic. 



Our interest in this letter lies in microstructures fabricated 
using electrostatic gates or etching that pattern a two dimen- 
sional electron gas in a semiconductor heterostructure, for 
example, GaAs/AlGaAs. Typically, the electronic transport 
mean free path is significantly larger than the dimensions of 
the device, and the electrons essentially travel ballistically 
across the microstructure. Their motion is governed by the 
shape of a smooth, self-consistent, steep-walled, confining po- 
tential which is often conceptualized as a quantum billiard. 

For many physical properties, the simplifying assumption 
that the dots' underlying classical dynamics are fully chaotic 
(hard chaos) has provided a good description of the experi- 
mental data [jjj H It has been used to justify various 
hypotheses from the applicability of random matrix theory 
(RMT) and random plane wave modeling (RPW) to statistical 
assumptions applied within semiclassical mechanics [^, ^, 
Indeed, chaotic systems manifest a large variety of universal 
behaviors. Furthermore, chaotic quantum dots are often qual- 
itatively very similar to diffusive ones provided the Thouless 
energy Eth is defined as hvp/L, where vf is the Fermi ve- 
locity and L is a typical dimension of the dot (as opposed to 
TiD/L 2 with D the diffusion constant). Consequently, most 
techniques, developed much earlier to study disordered met- 
als (diagrammatic approaches [Q], nonlinear sigma model [^|) 
and applied to disordered quantum dots [ |To[ |l2[ ], are ap- 
plicable to ballistic quantum dots. 

Nevertheless, unlike billiards, there are no known smooth 
potentials which are truly, fully chaotic. Unless designed oth- 
erwise for a specific purpose (such as measuring the weak 
localization lineshape [|l3[]) an odd-shaped, smooth potential 
genetically exhibits soft chaos, i.e. significant contributions of 
both stable and unstable motion. The general assumption of 
hard chaos is unfounded. 

As opposed to a genuine belief that the electrons' dynamics 
are strongly chaotic, the implicit assumption is that for many 
properties the distinctions between soft and hard chaos are 
more subtle than spectacular. This has been shown explicitly, 
for instance, for the fluctuation properties of Coulomb block- 
ade (CB) peak heights [|l4|], or for their correlations Jl5|]. In 
such circumstances, using a chaotic model allows for simpler 
analytic derivations without drastically altering the results. 

Our purpose here is to demonstrate that even this weaker 



assumption may, in some cases, be problematic; and that for 
some properties or measurements a strong sensitivity to the 
nature of the dynamics arises. Because the distinction be- 
tween the chaotic and integrable limits becomes most appar- 
ent at long (infinite) times, confining the electrons as long as 
practical within the dot holds the promise of leading to signa- 
tures of the dynamics. We shall therefore consider [zero tem- 
perature] ground states properties of well isolated dots, and 
find that ground state energies (whose second differences are 
probed by CB peak spacing measurements), and spin polar- 
izations are markedly affected by the dynamics. 

We proceed as follows. First, we give a general discus- 
sion of why, in principle, chaotic dots should be rather "atyp- 
ical", at least as far as ground state properties are concerned. 
Secondly, we consider in more detail a particular Hamiltonian 
model, namely a time-reversal non-invariant coupled quartic 
oscillator system, and show more quantitatively the relevance 
of the underlying dynamics. 

In the following, we assume a Fermi-Landau liquid descrip- 
tion of the quantum dot [ |To[ , |ll[ [lrj ], and therefore that the 
ground state energy E gr [N] is the sum of three terms 



E gr [N] = £ TF + £ lp + £ r] 



(1) 



£ TF [N] = (eN) 2 /2C is an electrostatic energy, £i p [N] is the 
sum of the single particle energies (SPE) of N independent 
particles, and £ rl is a residual interaction term. Specifically 



(2) 



with 6{ the SPEs corresponding to an effective potential 
V c ff(r) (that can be determined self-consistently within the 
electrostatic-like approximation) and /? = or 1 is the occu- 
pation number of orbital i with spin a = ±1, a ff = N). 
Furthermore, denoting ipi the eigenstate associated to ej, 

& = \ E fiff J drdr ' l^ 1 ")! 2 ^(r - r') |^(r')f 
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is the direct-plus-exchange, first order perturbation contribu- 
tion in terms of the weak screened Coulomb interaction V sc - 
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For the experimentally relevant electronic densities, V^ c is a 
short range function, and we shall model it as 



^ S c = ~S(r - r') 



(3) 



with v the total density of states (including the spin degen- 
eracy factor g s — 2) and £ e [0, 1] a parameter that can be 
related to the density of electrons (i.e. the parameter r s of the 
electron gas) and is in the range [0.5, 0.8] for many exper- 
iments. Within this zero range approximation, the residual 
interaction contributions read 
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with A the (local) mean SPE spacing and 



Mi. 



A I dr|^(r)| 2 |^-(r)| 2 



(4) 



(5) 



A is the dot area. The My are dimensionless and semiclassi- 
cal reasoning implies that their mean (for i ^ j) is unity. 

Note the crucial point that in Eq. (Q), only electrons with 
opposite spins actually interact. Aligning two spins decreases 
the residual interaction by a quantity of 0(A). There is thus 
a competition between the SPE term, which favors the occu- 
pation of the lowest orbitals, and the residual interaction term, 
which tends to align the electron spins. The relative strength 
of these two effects is governed by the dimensionless parame- 
ter £. If £ > 2, the ground state is completely polarized. This 
is the well known Stoner instability jl7[], which for instance 
is responsible for the ferromagnetic character of cobalt. 

Here £ is just less than one. Although full polarization is 
excluded, the proximity of the Stoner instability makes the 
ground state spin very sensitive to the fluctuations of the My 
and ej, which affects the fluctuation properties of the ground 
state [0, [Io[ [Tl], |l8| ], This is sometimes referred to as the 
"mesoscopic Stoner fluctuation". 

Now consider the diagonal term 



(6) 



M U =A / dr|V*(r)| 



Ma is the inverse participation ratio in position representation 
of the state ipi, which measures its extent of localization. Hard 
chaotic systems possess eigenfunctions that are the least local- 
ized in the sense that their Wigner transforms uniformly cover 
the energy surface Jl9[]. However, mixed systems are known 
to display various forms of phase space localization [p0||. This 
arises by various mechanisms. The most familiar one is asso- 
ciated with the quantization of invariant tori. It affects only a 
minority of states, but localizes them very strongly. Another 
mechanism discussed in [ p0[ ] is associated with the presence 
of partial transport barriers in phase space, the presence of 
which should be quite typical in mixed dynamical systems. 
Such partial barriers, if they are effective, should affect almost 
all eigenstates, but produce a lesser degree of localization. 



Similarly, it can be shown that the mean value of the off- 
diagonal terms Mij^i should be independent of the dynam- 
ics, and their fluctuations are extremely small for chaotic sys- 
tems (vanishing in the semiclassical limit [Q]). However, they 
would be of O(A) if significant phase space localization is 
present. We therefore see that chaotic systems, rather than be- 
having typically, are rather the limiting class of systems for 
which the interaction is the least effective. 

To explore the affect of soft chaos on the ground state prop- 
erties, we introduce a specific model. Let the effective Hamil- 
tonian H c ^ resulting from the lowest order, electrostatic-like, 
self-consistent calculation be (r = (x,y), r = |r|)) 

2 
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(7) 

Weyl operator ordering is assumed for the quantized version. 
The parameter b = 7r/4 is introduced so that the system has 
the symmetry of the rectangle instead of the square, a (A) is 
a convenient scaling factor chosen so that the mean number 
of states with energy smaller than E is given by N(E) = 
E 3 / 2 , regardless of the choice of A or k. A is the coupling 
between the oscillators. Finally, the parameter k breaks time 
reversal invariance (TRI). Note that for TRI systems, higher 
order terms in the ground state energy expansion, in particular 
the Cooper series, should be taken into account. Therefore, 
k is chosen such that TRI is completely broken. The specific 
form of the TRI breaking term has been chosen not to pertain 
to any particular physical realization, but to insure that the 
phase space portrait of the dynamics does not depend on the 
energy. By choosing A and k appropriately, various regimes of 
dynamics can be studied. In particular, we consider (A, k) — 
(+0.20, 1.00) [nearly integrable], (-0.20, 1.00) [mixed], and 
(—0.80, 1.00) [mostly chaotic]. Due to the TRI breaking term 
somewhat stabilizing the dynamics, the motion in the chaotic 
case is still not quite fully chaotic. 

The reflection symmetries lead to four irreducible represen- 
tations, which can be thought of as independent quantum dots 
with the same dynamics. We consider the four systems as 
an ensemble, which allows us to decrease our statistical 'er- 
ror bars' in the calculations. In fact, we increase the size of 
the ensemble even more by allowing A to vary ±0.02, which 
is enough to get nearly independent quantum eigenproperties, 
but small enough that the structure of the dynamics is essen- 
tially unchanged. Each statistical measure calculated within a 
given dynamical regime is thus the result of averaging over an 
ensemble of 12 similar quantum dots. 

For each parameter set, it is possible to compute for each 
symmetry class (e x = ±l,e y = ±1) the eigenvalues ti and 
eigenvectors ipi, from which the residual interaction terms 
Eq. (||) can be deduced (see [|1J for the numerical details). We 
first consider a few interesting statistical properties of these 
quantities, and afterward see how the spin distribution and 
ground state energy fluctuations are affected. 

In Fig. [j], values of sets of diagonal terms Ma are repre- 
sented for the symmetry class (+, +) in the various dynami- 
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FIG. 1 : Inverse participation ratio as a function of the orbital index 
for (+,+). From top down: (A,k) = (+0.20,1.00) [nearly inte- 
grable], (-0.20, 1.00) [mixed], and (-0.80, 1.00) [mostly chaotic]. 



cal regimes. The presence of very localized states in the nearly 
integrable and mixed regimes is immediately apparent. Curi- 
ously enough, it turns out that for all dynamical regimes, the 
IPR Ma is not very far from two for most states, but the distri- 
bution has a very long tail if the dynamics are not sufficiently 
chaotic. Consequently, it is reasonable to assume that the lo- 
calization mechanism at work here is associated with stable 
periodic orbits or tori. A more detailed analysis of the classi- 
cal dynamics reveals that partial barriers are actually present 
in (A, k) = (—0.20, 1.00), but that their position with respect 
to the symmetry lines of the system make them ineffective in 
changing the statistics of the Ma. Furthermore, and as can be 
seen in Table Q, the off-diagonal terms My (i ^ j) are also 
affected by eigenstate localization as their variance is signifi- 
cantly larger if the state i is localized. 

The difference between the distribution of My observed for 
the quartic oscillator system and that predicted for a chaotic 
system is essentially that a non-negligible number of very lo- 
calized states ipi give much larger diagonal terms, and much 
larger variation of the corresponding off-diagonal terms. The 
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FIG. 2: Integrated peak spacing distribution with ( = 0.8 for 
RMT/RPW prediction [Q| (dotted line), mostly chaotic (solid), mixed 
(dash-dot) and nearly integrable (dash) regimes. 



question is how does this affect the ground state properties. 
The answer for given (A, k) and number of electrons in the 
dot N follows from the &j and M«. To compute the en- 



ergy E[{f?}], we use Eqs. (|lj,^|) for the various occupa- 
tions {/f} such that ^2 icr f[ = N. The ground state fol- 
lows by selecting the occupation sequence minimizing this 
energy. Varying N in the range [100, 200] for each param- 
eter set, one constructs a distribution of total spins, occu- 
pancies, and second differences in the ground state energy 
£[N] = E SI [N + 1] + E gI [N - 1] - 2E gI [N]. £[N] is exper- 
imentally accessible by measuring CB peak spacings, and we 
therefore refer to it as a "peak spacing". The various dynami- 
cal regimes' peak spacing distributions are shown in Fig. ||. 

We see that the observed distributions are strikingly differ- 
ent from the RMT/RPW predictions (shown on the same plot) 
that should apply for fully chaotic systems. Even our most 
chaotic case shows significant deviations, and this increases 
considerably as the system moves away from hard chaos. The 
peak spacings fluctuate much more, and very large spacings 
appear for the mixed and nearly integrable regimes. More- 
over, as shown in Table |j[ larger spins become significantly 
more probable. 

We next ask how is such a drastic change made possible by 
the presence of a moderate number of very localized states. 
Let's consider one specific example in detail. Fig. || shows the 
succession of orbital fillings for a range of N. Here, two of the 
single particle states(i = 64, 66), are highly localized. Since 
the system is not chaotic, there is also less level repulsion, 
which allows for the clustering of levels around £64 or egg that 
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Mu/gs > I snp 


0.097 


0.108 


0.070 


Mu/gs < 7 inf 


0.024 


0.023 


0.009 



TABLE I; Conditional variance of the interaction terms My with 
< \i — j\ < 10, i > 51. Rows top down: i) dynamical case, ii) su- 
perior limit, iii) inferior limit, iv) conditional variance with localized 
orbitals, and v) conditional variance with delocalized orbitals. 
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+0.20 


-0.20 


-0.80 


RMT/RPW 


P( S = 2) 


0.13 


0.16 


0.07 


0.01 


P(s = 5/2) 


0.08 


0.10 


0.02 


0.00 


(8s) 


0.51 


0.54 


0.38 


0.23 



TABLE II: Probabilities P(s = 2) , P(s = 5/2) to find a spin two 
(even N) or five halves (odd N) ground state, and average value (Ss) 
of the ground state spin augmentation (Ss = s or (s — 1 /2) for even 
or odd number of particles, respectively), for the various dynamical 
regimes (values of A) with k — 1.0 and £ = 0.8. The last column is 
the RMT/RPW prediction 
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FIG. 3: Successive filling of the (+, +) orbitals i — 63 to 69 as 
the number of particles in the dot goes from TV = 129 to 135 for 
(A, re) = (+0.20,1.00) and < = 0.8. The spacing between the 
horizontal lines are proportional to the actual level spacings. The 
numbers in the right column are the corresponding IPRs. 



space localization in their single particle properties, whereas 
this point is debatable with respect to the dots of [fit]. How- 
ever, it may be more interesting to address this question from 
the opposite point of view. Since creating dots away from the 
hard chaos limit leads to behaviors that are qualitatively differ- 
ent from those predicted using chaotic or diffusive modeling, 
if the interest is in devising a dot to perform some particular 
function, such as spin manipulation, for instance, it is in the 
soft chaos regime that richer behavior involving large fluctua- 
tions of ground state energies and spins will be found. 

We appreciate discussions with H. Baranger and G. Usaj. 
This work was supported in part by NSF Grants # PHY- 
0098027 and DMR-0 1 03003 . 



would otherwise be improbable. We see that none of the or- 
bital occupations follow the simple "up/down" scenario char- 
acteristic of non-interacting systems. As would seem natural, 
the very localized orbitals i = 64, 66 remain singly occupied 
across many values of N. There is an additional non-intuitive 
feature, namely that the other orbitals also prefer single occu- 
pancy despite not being particularly localized. 

This behavior derives from the following mechanism. Sup- 
pose for whatever reason, a given orbital i* is singly occupied, 
say with spin +. The M«* have to be included in the total 
energy for all orbitals i if (i, — ) is occupied, but not if only 
(i, +) is. This can make the energy cost of occupying an or- 
bital with spin down (£A/g s ) higher than doing so with spin 
up. Assuming Ni such singly occupied orbitals already exist, 
and adding the term ((A/ g s )Mu associated with the interac- 
tion between the (£, +) and the (£, — ) particles, on average it 
turns out that (even for a non-localized state) the residual in- 
teraction energy cost of doubly occupying some orbital i is 
((JVj + 2)(A/g s ). For typical values of £, this is larger than a 
mean spacing as soon as TV; > 1. Consequently, as illustrated 
by Fig ||, the localized orbitals will not only remain singly oc- 
cupied, but also have a tendency to polarize the electrons in 
the other nearby orbitals. Lack of level repulsion and larger 
fluctuations of the Mij will further enhance such effects. 

To conclude, we have shown that for the non-TRI quar- 
tic oscillators the eigenfunction statistics behave differently 
than those from a RMT/RPW approach. The significance is 
correlated with the degree of chaos (or lack thereof) in the 
underlying classical dynamics. Due to the proximity of the 
Stoner instability, strong effects arise in ground state spin po- 
larizations, occupancies, and energies of the corresponding 
"model" quantum dot. The quartic oscillators more fairly rep- 
resent a generic experimental dot than the hard chaos assump- 
tion. The hard chaos assumption leads to predictions that are 
qualitatively and quantitatively incorrect. 

Finally, such considerations should affect the understand- 
ing of realized dots, which needs to be discussed on a case- 
by-case basis. For instance, the dots used in [ p2| ] are certainly 
far from chaotic. They should show a large degree of phase 
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